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Abstract: We calculate the rate of energy loss from compact astrophysical objects due to 
a scalar field in screened modified gravity models of the chameleon, dilaton and symmetron 
types. For binary systems, this is related to the time drift of the orbital period comple- 
menting the effect due to gravitational wave emission in General Relativity. The scalar 
energy flow is sensitive to the time evolution of the scalar field on cosmological scales and 
the screened gravitational charge of the massive bodies, eventually leading to a dependence 
on the range of the scalar interaction on cosmological scales in the present Universe. Using 
the Hulse- Taylor binary pulsar, the double pulsar PSR J0737-3039 and the pulsar-white 
dwarf system PSR J1738+033, we find the most stringent observational bound on the range 
A of the scalar field A < 1 Mpc. 
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1. Introduction 

The acceleration of the Universe [1,2] may be due to the existence of modified gravity on 
large scales [3]. In most known scenarios, this possibility is realised by augmenting the 
gravitational sector with a scalar field. Scalar-tensor theories lead to detectable effects 
in the solar system which have been constrained drastically by the Cassini probe [4]. In 
practice, this requires that the scalar field must be screened in dense environments. Such 
models abound and fall within three distinct categories which can be easily understood 
using the Lagrangian 

C D -^W) 2 - + , (1.1) 

2 2 mpi 

corresponding to small variations 5(f> of the scalar field around the background value <j)Q. 
The scalar field couples to the trace of the matter energy momentum tensor 5T. The 
scalar force is screened by the Vainshtein mechanism [5] when Z^^) is large enough 
that the coupling of the normalised field f3(4>o)/Z 1 / 2 ((j)o) is small enough. The chameleon 
mechanism [6-8] occurs when the mass m(0o) is large enough to suppress the range of 
the scalar force. Finally, the symmetron or dilaton screenings are such that (3(4>o) itself is 
small [9-12]. In all these cases, the background value (fto depends on the environment and 
the screening mechanisms occur in the presence of dense matter. 

The screened models of modified gravity of the chameleon, dilaton and symmetron 
types can be described in a unified way [13] using a tomographic method whereby the 
potential and the coupling function to matter as functions of the scalar field can be com- 
pletely reconstructed from the time dependence of the minimum of the effective potential 
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in the cosmological background since Big Bang Nucleosynthesis (BBN). This allows one to 
parameterise screened models using the value of the scalar mass now rrioo and its coupling 
strength ftoo. In particular, imposing that solar system tests of gravity are evaded and that 
the Milky Way is screened to avoid too much disruption of the dynamics of its galactic 
halo, one deduces that the scalar mass must be such that > 10 3 Hq [14,15]. This bound 
is crucial as it gives the range of the scalar interaction on cosmological scales and therefore 
restricts the possible astrophysical effects of screened modified gravity. In particular, this 
implies that future galaxy surveys will not probe screened modified gravity on scales larger 
than 10 Mpc. In this paper, we tighten up this bound using the time drift of the orbital 
period P of binary systems. These measurements place one of the tightest constraints on 
General Relativity (GR) and provide indirect evidence for the existence of gravitational 
waves. 

Here we calculate the scalar energy change in a finite volume per unit time for screened 
modified gravity models with canonical kinetic terms, therefore excluding k-essence and 
Galileon type models subject to the Vainshtein mechanism [16, 17], and find that it depends 
on both the gradient of the scalar field and its time derivative. The gradient is generated 
by the binary system and leads to a 1/r decrease of the field in the far away regime. 
As the gravitational charge of screened objects such as neutron stars and white dwarfs 
is tiny, the usual monopole, dipole and quadrupole contributions to the scalar radiation 
are negligible. The contribution due to the time dependence of Newton's constant is also 
irrelevant. The only time dependence which is of relevance to the energy loss is entirely due 
to the cosmological evolution of the field background value. In the absence of cosmological 
dynamics, hardly any scalar energy flow would be induced as the screened gravitational 
charge of each object is minute. It turns out, and this is the main result of this paper, that 
the induced P is small enough when the range of the scalar force is smaller than 1 Mpc, . 

The paper is arranged as follows. In the first section, we evaluate the rate of change of 
the scalar energy in a classical setting. We then reassess the argument for binary systems in 
the relativistic case. Using tomography we relate the change of energy to the cosmological 
range of the scalar force. Finally we apply our results to the cases of binary pulsars and 
pulsar-white dwarf systems and deduce a bound on the range of the scalar interaction on 
cosmological scales. 

2. Classical Scalar Flow 

As a first step, let us analyse the radiation of scalar energy in scalar-tensor models where 
the background metric is assumed to be Minkowskian. We are interested in a scalar-tensor 
theory where the scalar field is screened in the presence of dense matter. This theory is 
characterised by its bare potential V(cf>) and the coupling function A(cf>) 

S = J d'x^i^- - - V(cf>)) + S m (A,A 2 (<f>) 9lu/ ) (2.1) 

where the particles represented by the fields ipi interact with the rescaled metric g^v = 
A 2 (<j))g liv . As the kinetic terms are canonically normalised, we exclude from our analysis 
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all models subject to the Vainshtein mechanism which require a special treatment [16,17]. 
One of the salient features of the models described here is that the effective potential of 
the scalar field in the presence of matter is modified and becomes 

V cS (4>) = V(<f>) - (A(cf>) - 1)T (2.2) 

where T = — p + 3p is the trace of the energy momentum tensor of matter. Notice that we 
have subtracted the contribution from matter, hence the presence of the T term, in order 
to separate the scalar energy from the matter energy. In this section, we concentrate on 
the pressureless case where p is time-independent and localised within a finite volume. We 
focus on the energy inside this domain which is simply 

E = j^x^ + ^ + V^))- (2-3) 

In this approximation, the scalar field obeys the Klein Gordon equation which is simply 
given by 

Multiplying by (j> and integrating over a domain D with boundary dD sufficiently far from 
the compact objects in D, we get 



dtJ D dX[ 2 



2 + V eS ) = I d^Ac/). (2.5) 

D 



Integrating by parts 



f d 3 x4>A(f>= [ d 3 x(di(<j)di0) - di<j)di(j)) (2.6) 
Jd Jd 

from which we get the energy loss due to scalar radiation 

d 2 x4>di4>m (2.7) 



dE 
~dt 



dD 



where we have used Stokes' theorem and n; is the outward pointing normal to the boundary 
(n 2 = 1). The flux term comes from the current 

Ji = 4>di<p. (2.8) 

In the following, we will see that the energy flow is universal for screened objects in modified 
gravity and depends essentially on the time drift of the scalar field on cosmological distances 
due to the time dependence of the Cold Dark Matter (CDM) energy density. 



3. Energy flux from Binary Systems 

The previous results cannot be directly applied to astrophysical situations of interest in 
as much as the relativistic effects present in binary systems containing neutron stars are 
not taken into account. One must also consider the fact that the background metric is 
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not Minkowskian and converges to an accelerating Friedmann-Roberston- Walker (FRW) 
Universe at infinity. In practice, we are interested in time scales corresponding to the orbits 
of binary systems, i.e. a fraction of a day, which are much smaller than the age of the 
Universe. Hence we can safely approximate the asymptotic FRW Universe with Minkowski 
space-time when studying the scalar radiation from binary systems. In the relativistic case, 
a proper formulation of the conservation of gravitational and matter energy is more difficult 
to obtain. Here we will make use of the gravitational energy momentum pseudo-tensor t^ u 
defined by Landau and Lifschitz [18] whereby the total conservation law reads 

d,((-g)(T^ + T^ + t^)) = 0. (3.1) 

The matter energy momentum tensor is while the one of the scalar field is T£ v . This 
allows one to identify the total energy density as 

T ° = g(T° m + + tg) (3.2) 

and the energy flow vector far a way from the binary system as 

n = %+ti (3.3) 

where asymptotically we have g = 1 as the Universe converges to Minkowski space-time. 
Integrating over a large volume surrounding the matter distribution and comprising the 
astrophysical objects which are radiating, we find that the time dependence over time scales 
much shorter than the age of the Universe (hence we can consider the scale factor a® = 1 
now and during the whole radiative process) of the total energy is given by 

^ = F <t> + F giav (3.4) 

where the total energy inside the volume D surrounding the binary system is identified 
with 

E= d 3 xg(T° m + T^ + t° ) (3.5) 

JD 

and the energy fluxes emanating from the binary system can be separated into two com- 
ponents 

IdD ' JdD 
The gravitational flux corresponds to the gravitational wave emission while the scalar flux 
results from the dynamics of the scalar field. Moreover we have of course 



/ T^n u F grav = / tlm (3.6) 

JdD JdD 



t; = r = m. (3.7) 

In the following, we will concentrate on the flux due to the scalar field and compare it 
to the flux from gravitational waves. Moreover we will always assume that the Yukawa 
suppression of the scalar field due to its finite mass moo far away from the binary system 
is negligible. This requires that the range of the scalar interaction A = is much 

larger than the typical distance between the stars of binary systems, i.e. a few kpc. In 
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practice, we will find that A must be less than 1 Mpc, complying with our hypothesis in 
the 1 kpc < A < 1 Mpc range. 

For relativistic stars in their interior, the pressure p cannot be neglected and T = 
—p + 3p [19-21]. When the scalar field is screened in a dense object, the effective potential 
inside has a minimum <p c which depends on the details of the stellar model through the 
equation of state p(p). On the other hand, very far away from the compact object, the 
effective potential is sensitive to the pressure-less density of Cold Dark Matter (CDM) and 
has a minimum (j)^. For neutron stars, the Newtonian potential at the surface of the star 
<3?jv = G^M/R, where M is its mass and R its radius, can be as large as 0.2, in natural 
units where 8ttGn = V m pi an ^ m Pi ~ 2 • 10 18 GeV is the reduced Planck mass, implying 
that relativistic effects must be taken into account inside the star. This potential is just 
below the instability threshold <3?tv = jg, for stars with a constant equation energy, where 
the trace contribution to the effective potential leads to an instability [19,20]. For values of 
$tv below this threshold, we can describe the far away field generated by a single screened 
object as 

4>{r) = 4>oo + Q m Pi— ^— (3-8) 
r 

where Q is the gravitational charge of the scalar field and we consider distance scales 
much lower than the inverse mass of the scalar field allowing one to neglect the Yukawa 
suppression term due to the mass of the scalar field at infinity 1 . This form of the solution 
follows from the fact that far away from the compact object the metric can be taken to be 
Minkowskian and the scalar field satisfies the Laplace equation Ac/) = 0. We are working 
in the quasi static approximation where the time derivatives are much smaller than the 
spatial gradients. A small time variation of the field is induced by the cosmological time 
dependence of the field value at infinity 0oo(*)- I n the Minkowskian approximation, the 
field generated by a single screened object [7, 13] is given by 

<f>(r) = <l>o (t) + (<f> c -<f> 00 (t))- (3.9) 

Here cp c is the value inside the dense body which is time independent and is obtained by 
minimising the effective potential taking into account both the matter density and the 
pressure inside the compact object. The gravitational charge Qa of each screened object 
A is not universal and depends on the Newtonian potential at its surface &a 

Qa = (3.10) 

The gravitational charge of screened objects Qa is independent of the bare charge at infinity 
/?oo = rapi 9 \ <f>=4, ao - The lunar ranging experiment places a bound on a violation of 
the equivalence principle in the earth-moon system implying that [7] 

Qe < ID" 7 (3.11) 



1 In the chameleon case [7], this was originally identified as 6/3^ where AR is the size of the thin 
shell. This result has been generalised to all screened modified gravity models of the non-Vainshtein type 
in [13,22-24]. 
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for the earth. For much denser objects, such as white dwarfs and neutron stars, the 
gravitational charge is even smaller Qa <C 1CP 7 . This confirms that white dwarfs and 
neutron stars are screened astrophysical objects. 

As long as the value of the field at infinity is considered to be constant, the scalar 
radiation in scalar-tensor theories has been extensively studied [25-28] . In the case of binary 
systems, the resulting energy flow, which follows from the motion of the astrophysical 
objects, can be either monopolar, dipolar or quadrupolar [25]. The case of monopolar and 
quadrupolar radiation is interesting as [25] 

(^) mono - quadm oc f 2 (Q A , Qb)^^ (3.12) 

dE 

where — jp 1 is the emission due to gravitational waves and /2(Qa,Qb) is bilinear in the 
charges of the two bodies Qa,b- As Qa,b *C 10~ 7 for neutron stars and white dwarfs 2 , 
the monopolar and quadrupolar fluxes are negligible compared to the gravitational wave 
emission. The most stringent test of scalar emission when the field at infinity has negligible 
time variation comes from the dipolar emission rate [29] for two bodies with different 
gravitational charges, e.g. for a binary system comprising a neutron star and a white 
dwarf [25] 

, ^dipolar 2tt M A M B 2 2ttG n (1 + Q a Qb)M 1 + e 2 /2 

( ^r } = ~i+q a Qb m 2 p (i _ e 2 )5/2 ( 3 - 13 ) 

where Ma,b are the masses of the objects, M = Ma + Mb, P the orbital period and e the 
eccentricity. Observations of white dwarf-neutron star systems such as PSR J1738+033 [29] 
imply that the amount of dipolar radiation is negligible when Qns *C Qwd *C 1CP 7 as 
is the case of neutron stars and white dwarfs. All in all, the monopolar, dipolar and 
quadrupolar amount of radiation from a screened system of binary stars is negligible. 

The rate of energy flow ^ can be compared to one of the salient characteristics of 
binary systems: the time drift of the orbital period It is commonly assumed that the 
energy flow due to gravitational waves and scalar field radiation has a direct effect on the 
binding energy of the binary systems (for a more thorough discussion, see [25]) which leads 
to a non-zero value for the time variation Indeed, using the fact that the period P 
satisfies Kleper's third law to leading order a 3 /P 2 = constant and that the gravitational 
energy scales as Eg ~ a -1 , we find that the energy flow due to the scalar field induces a 
change in the orbital period of a binary system as 



dlnP^ _ 3 dE 



(3.14) 



dt 2 E G dt 

In the following, we will transcribe our results on scalar radiation in terms of the time 
variation or the orbital period 



2 It is interesting to notice that screened models of modified gravity realise a weakening of the gravitational 
charge for dense objects Q <C fioa which is due to the existence of a minimum of the effective potential in the 
presence of matter. When the bare potential is negligible and the coupling function is of runaway type such 
as A((f>) — exp(—c(f> 2 ), the field 4> is driven to large values in dense environments where the gravitational 
charge becomes large, an anti-screening behaviour which was particularly emphasized in [26]. 
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Another effect affecting ^? arises when Newton's constant has a time dependence, as is 
the case if one considers the Einstein-Hilbert action with an explicit time dependent G/v (t) . 
This leads to a time variation of the period of binary systems due to the time change of 
the gravitational constant [30]. This contribution can play an important role and satisfies 
present observational bounds [31-33] on deviations from General Relativity provided [34] 

&|< 0.02# (3.15) 

Gn 

as evaluated in the recent Universe. In scalar tensor theories such as the ones leading to 
screened modifications of gravity, a time drift of Newton's constant is present in the Jordan 
frame where matter couples minimally to the metric (and where the orbits are therefore 
calculated), it reads 

1^1 = 2/3^ (3.16) 

where (j)^ is evaluated in the Einstein frame. We will see that the resulting bound is easily 
satisfied as soon as the screening of the Milky Way is imposed on the model parameters. 

On the other hand, a new effect involving the time drift of cf)^ can lead to a non- 
negligible amount of scalar radiation, resulting from the 1/r profile of the scalar field on 
large scales 3 . Let us now consider a binary system where the typical distance between the 
stars is much larger than the radii <C D. The study of the coupled dynamics involving 
the scalar field and the metric in the relativistic regime requires numerical methods. In this 
paper, we will use a simplified approach whereby the pressure inside the star is taken into 
account but relativistic corrections to the metric are neglected. This can only be considered 
as a first approximation to the dynamics, valid at the ten percent level, which will be enough 
to capture the salient features of the scalar flow from compact objects. More precise results 
would necessitate the study of the time evolution of the compact system taking into account 
the non-trivial metric and would have to be carried out numerically. This is beyond the 
scope of this work. A reasonable approximation for the very far away solution from the 
two objects is given by the superposition of the two screened solutions obtained for a single 
object 

<f>(r,t) «0i(r- rx(t)) + <h(r- r 2 (t)) - (3-17) 

where J"i 2(t) are the time dependent positions of the two bodies. This is valid in the case of 
non-relativistic motions where |ri 2I <C 1 allowing one to neglect time derivatives compared 
to spatial gradients. In the far away approximation we have 

4>(r,t) « ^(t) - <^(t)^I±^ (3.18) 

where we have assumed that 4> c <C cj)^ as the density at infinity is very small compared 
to the density inside the compact objects. This implies that, to leading order in 1/r, the 
current is given by 

Ji = 4>oo<Poo^-^2~^ n i- (3.19) 



3 For scales smaller than the range A 



-7- 



Notice that the current is non-vanishing only because of the time variation of the value at 
infinity. It is a feature of screened modified gravity models that local properties, such as 
deviations from Newton's law and here the energy flow, are sensitive to the behaviour of 
the field at infinity, a Machian property. This induces a finite contribution to the energy 
flow due entirely to the scalar field and which is given by 

dE 

= 4^00^00(^1 + R 2 ), (3.20) 

depending on the dynamics at infinity only. 

The time drift of the energy E comprises two terms: the gravitational contribution and 
the scalar field one. The first one is due to the emission of gravitational waves. The second 
one is new, in particular here where its main component results from the time variation of 
the scalar field at infinity. As customary in the literature, we shall express the scalar flux 
as a contribution to the change of the orbital period of the binary system which becomes 

^ = 6^00^00(^1 + ^2)^ (3.21) 

where the gravitational energy is approximated by 

M X M 2 

E G = ~ Te , 2 ■ (3-22) 

Of course, the contribution we have calculated is only a new physics contribution to P 
and its value can only be compatible with experiments if it is smaller than the discrepancy 
between the General Relativity (GR) contribution due to gravitational waves and the actual 
measurements: 

1 dt 1 - 1 dt dt '' [ ' ' 

For the case of binary systems such as binary pulsars or pulsar-white dwarf systems, this 
new physics contribution must be tiny implying strong constraints on the modified gravity 
models. 



4. Modified Gravity Tomography 

The time evolution of the cosmological vacuum value 0oo for all screened modified gravity 
models is known and is given by [14] 

9Ax,fim0^f-#ompi (4-1) 



at m 



■DC 



where moo is the mass of the scalar at infinity and /3oo the coupling at infinity. It involves 
the Hubble rate now and the matter fraction f2 m o- We can also have access to the value of 
the field at infinity 0oo thanks to the tomographic property of these models. Indeed, the 
dynamics of screened modified gravity models are entirely determined by the time evolution 
of the minimum of the effective potential in the course of the cosmological history. The 
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value of the field in a dense and pressure-less medium whose density corresponds to the 
CDM matter density p for a value of the scale factor a is given by [14] 

~ <Pc 3 f a (3(a) 
mpi mpi J a . n . am 2 (a) 

where 4> c is the minimum value in a dense object whose density corresponds to a small value 
of the scale factor a m \. The time evolution /3(a) and m(a) defines the particular model of 
modified gravity under study. As already stated, in practice <\> c <C 4>(a) as soon asa> a^. 
In particular, this implies that the value of the field in the cosmological background now 
is given by 

H 2 

</>oo = 9fW/3oc^ (4.3) 

where the cosmological evolution of the coupling and the mass are m(a) = m 00 f(a) } /3(a) = 
Poo9{a) and I = da a Ap^ is a coefficient of order one which is model dependent. We 
have taken a- m \ ~ as the precise value is irrelevant as long as the integral / converges. 
From this we get that 

This depends on the ratio tti^/Hq in a crucial manner. From the comparison between the 
experimental results and the GR calculation of the gravitational wave emission, a bound 
on l^^l can be deduced. This in turn leads to a lower bound on and an upper bound 
on the range A = of the scalar force on cosmological scales. 

Let us give an example of this reconstruction procedure (more examples can be found 
in [13]). An interesting class of models corresponds to f(R) gravity in the large curvature 
regime where 

/3(a) = p oo = -L (4.5) 

and 

m(a)=m 00 ( 4 "^ + "y " 3 )(^)/ 2 (4.6) 
4J I AO + s 'mO 

where Q\ ~ 0.7 and fi m o ~ 0.3. Here n is an index which is unspecified. More complex 
examples like symmetrons or dilatons can be easily considered too. 

This mass bound deduced from a bound on \^f-\ can be compared to the one obtained 
by imposing that a galaxy such as the Milky Way must be screened. This guarantees that 
the field inside the galaxy is constant and that solar system bounds on modified gravity 
are satisfied. This also prevents strong disruption effects on the dynamics of the galactic 
halo although the extent of this phenomenon would have to be ascertained by numerical 
simulations. The screening condition in the Milky Way can be formulated as 

I<Ag-<Aoo| <2/3oo<I> G (4.7) 

where the Milky Way has a Newtonian potential at its surface equal to <3>g and a value at 
the minimum equal to <Pg- This can be reexpressed as 

ln m0 ^- Hda-^L-KQo. (4.8) 
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In the case of the Milky Way ac = 10 -2 corresponding to a density inside the galaxy which 
is 10 6 larger than the cosmological density and ao = 1 assuming that the Milky Way is 
surrounded by the cosmological vacuum. As long as the integral fjf° da a fj = 0(1), and 
upon using that the galactic Newtonian potential is &g ~ 10 -6 , we get that [14,15] 

> 10 3 . (4.9) 



In the following, we will see that a more precise bound can be obtained from the timing of 
binary systems. 

We are now in position to verify that the contribution to the change of the orbital 
period due to the time variation of Newton's constant is negligible. Indeed for (3^ = 0(1) 
and upon using fl4.1|), we have 



= 9/?^Q m0 ^^o < 10- 6 i?o (4.10) 



implying that the bound ( |3,15 ) is easily satisfied. We now turn to the new contribution to 



scalar field radiation (M.4h. 



5. Phenomenology 

The time variation of the period P that we have obtained can be applied to some very 
precise observations of binary systems. This will allow us to deduce a bound on the mass 
moo of the scalar field responsible for screened modified gravity. Moreover this constraint 
is competitive with the one deduced from the screening of the Milky Way in as much as 
it does not rely on the very complex dynamics of galaxies. It only depends on the better 
understood physics of binary systems and can be taken as a stronger constraint than the 
one deduced from the screening of the Milky Way. We will consider three cases where 
P exp has been measured and compared to -Pgr- The contribution from the scalar field 
must be compatible with these measurements and therefore be bounded by the uncertainty 
between the experimental value and the GR result. In all cases, the input parameters are 
the semi-major axis a, the masses Ml 2 and the radii R± 2. In the neutron star case, the 
radii are not known precisely but reasonably inferred assuming that the equation of state 
of nuclear matter is the one of an ideal gas, this is a crude assumption which suffices for 
the type of bound we are looking for. Moreover the masses Mi 2 of the two stars can 
be inferred from the knowledge of two observables: the rate of advance of periastron Co 
and the timing parameter 7. It turns out that the calculated values of these observables 
is independent of (f)^ and deviate from their GR values in a way which depends on the 
gravitational charges Qi t 2 quadratically. For charges Qx$ <C 10~ 7 as obtained for neutron 
stars and white dwarfs, the deviation from GR is negligible and one can therefore infer the 
two masses Mi i2 and use their values deduced in GR [27]. 

The loosest bound is obtained with the double Pulsar PSR J0737-3039 [35] where 
a = 878959796.549 m, Mi = 1.3381(7)M , M 2 = 1.2489(7)M Q , P = 0.10225156248(5) day, 
and |P cxp | = 1.252(17)10- 12 at the la level. The GR value is |P GR | = 1.24787(13)10^ 12 . 
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The difference is at the 5 1CT 15 level. We will take R 1 ~ 13.72 km and R 2 ~ 14.04 km. As 
a result we obtain that 

^Idouble pulsar > 109 ({U) 
-"0 

where /?oo and / are model dependent. This implies that structure formation is affected on 
scales of order 

k > 0.26 / 1/4 /4 /2 Mpc" 1 (5.2) 
or in real space the range of the fifth force is 4 

Adouble pulsar < 38 /' 1/4 /3" 1/2 M P C. (5.3) 

A better bound can be obtained using the Hulse Taylor pulsar PSR B1913+16 [36] consist- 
ing of a pulsar and a companion neutron star. In this case we have a = 1950100 km, M\ = 
1.4398(2)M , M 2 = 1.3886(2)M , P = 7.751939106 hours, and |P cxp | = 2.425(1)10^12 at 
the la level. The GR value is |P GR | = 2.40458(4)10~ 12 . The difference is at the 2 10~ 15 
level. We will take R\ ~ 13.39 km and R 2 ~ 13.55 km. Hence we obtain that 

^|Hulse-Ta y lor> 257 / 1 / 4 /3^ 2 . (5.4) 

This implies that structure formation is affected on scales of order 

k > 0.06 / 1/4 /4 /2 Mpc- 1 (5.5) 
or in real space the range of the fifth force is 

AHulsc-Taylor < 16 I' 1 ' 4 ^ 2 M P C. (5.6) 

The most stringent bound can be obtained with the measurements of the pulsar- white dwarf 
binary system PSR J1738+033 [29]. We have a = 1738982401 m, M P = 1.46 + 0.06/ - 
O.O5M , M WD = 0.181(+0.008 - 0.007)M Q , P = 0.3547907398724(13) day, and |P cxp | = 
25.9 ± 3.210~ 15 at the la level. The GR value is |P G r| = 27.7t}JlO- 15 . The difference is 
at the 2 10~ 15 level. We take R P ~ 13.33 km and R WD = 0.037(+0.004 - 0.003)/?© km. 
As a result we obtain that 



rJ "[pulsar-white dwarf > 3623 I 1 ^^ 2 . (5.7) 

Mo 

This implies that structure formation is affected on scales of order 

k > 0.87 I 1/4 PU 2 Mpc- 1 (5.8) 
or in real space the range of the fifth force is 

Apulsar-whitc dwarf < 1-14 r 1/4 ^ 1/2 M P C. (5.9) 

Hence we expect that modified gravity effects on large scale structure can only be felt at 
the Mpc level, in the non-linear regime of structure formation which can only be probed 
using N-body simulations. 



4 The scalar interaction is Yukawa suppressed on distances larger than the inverse mass of the scalar field 
on very large scales A = m^} . 
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6. Conclusion 



We have analysed the time drift of the orbital period of binary systems due to the flow of 
a scalar field modifying gravity on large scales. The rate of change of the orbital period 
due to the scalar field is related to the cosmological time evolution of the scalar field and is 
therefore linked to the range of the scalar force on cosmological scales. We have found that 
the observation of the pulsar- white dwarf system PSR J1738+033 places the strongest 
constraints on the scalar range which has been found to be of the order of 1 Mpc. This 
implies that observations of modified gravity effects on astrophysical scales can only occur 
in the sub-Mpc regime where numerical simulations taking into account both modified 
gravity and baryons will have to be envisaged [37]. This is left for future work. 
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